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Based on standard k-p (8 x 8) multiband Hamiltonian, we have deduced an explicit analytical 
expression for the Rashba-coupling parameter which clarifies its anomalous behavior for heavy holes 
{hh), gated in quasi-two-dimensional (Q2D) systems, by letting grow the density. Our modelling 
remarkable better agrees with experimental results in comparison with earlier theoretical models, 
while recovers the expected cubic dependence on the quasi-momentum. For quasi-one-dimensional 
(QID) hh systems, we have formally derived an effective Rashba Hamiltonian with two competitive 
terms on the quasi-momentum, a linear term and a cubic one as predicted from suitable approxima¬ 
tions to the Q2D scope. The Rashba-coupling parameters also behave anomalously and qualitatively 
support recent experiments in core/shell nanowires. Furthermore, they exhibit an essential asymp¬ 
totic discontinuity in the low density regime as a function of the lateral confinement length. For hh, 
we present closed schemes to accurately quote the Rashba-coupling parameters both for the Q2D 
and QID systems, which become unprecedented for holes. 

PACS numbers: 71.70.Ej, 73.21.Hb, 73.21.Fg, 85.75.-d 


One of the biggest challenges in Spintronics [l[ is to 
manipulate efficiently spin currents in semiconductors 
systems without external m^netic fields. Rashba spin- 
orbit interaction (SOI-R) [IjU is amoM the most promis¬ 
ing mechanisms for achieve this [J-l7|. SOTR for the 
heavy holes (hh) case is different from that for electrons 
or light holes (Ih). In quasi-two-dimensional (Q2D) sys¬ 
tems, observed in experiments since 1984 Q, the SOI- 
R coupling Hamiltonian exhibits an atypical cubic de¬ 
pendence in the quasi-momentum (k^) obtained by Ger- 
chikov and Subashiev [^, and Winkler et al. [IMl us¬ 
ing group theory arguments. 

The long data assumption that the induced SOTR spin 
splitting in Q2D systems at zero magnetic field rises with 
electric field, responsible for the inversion asymmetry 
of the confining potential, firmly validated for electrons 
, was definitel y a bated when the opposite was demon¬ 
strated for holes [l0|. Importantly, it has been reported, 
both numerically and experimentally, the SOTR coupling 
parameter anomalous decreasing with a non-zero electric 
field for Q2D hh in accumulation-layer-like single het- 
13, HI- Recently, Habib et al. [h 


erostructures 


have 


shown interesting measurements, with a similar behav¬ 
ior at fixed density, but tuning the external electric field 
perpendicular to the Q2D gas. Notwithstanding these 
achievements, several features in the anomalous SOTR 
coupling effect remains cumbersome and inspired us to 
address a complementary study. 

On the other hand, despite the efforts focused to study 
SOTR in quasi-one-dimensional (QID) holes systems, 
due to their appealing applications and phenomenology 
IMl, some relevant topics remains incomplete still or 


deserve further attention. SOTR model for QlD-hh sys¬ 
tems suggested by Governale and Ziilicke [15|, and in¬ 
tuitively based on the behavior of the Dresselhaus term 
for Q2D electrons [i^, leads to a Hamiltonian linear in 
the quasi-momentum. However, Zhang and Xia [l^ sug¬ 
gested that additional competitive high-order terms may 
be needed to explain their finding regarding the spin 
splitting energy of hh in cylindrical quantum wires un¬ 
der external electric fields. Tikewise, from our calculation 
of hh dispersion laws in QID systems 0: patterned in 
Q2D hole gases by repulsive bias, we demonstrated the 
existence of a cubic-in-the-quasi-momentum term. Re¬ 
cently, Chesi et al. [l^ proposed a clever model, replac¬ 
ing the quasi-momentum operator along a constricted 
direction of a Q2D-hh system by its expectation value. 
They obtained two competitive terms, a linear and a cu¬ 
bic one, that explain the anomalous sign of the spin po¬ 
larization filtered by a quantum point contact (QPC), 
as observed in magnetic focusing experiments, as well as 
the crossing or anticrossing of spin-split levels depend¬ 
ing on subband index and the magnetic field direction 
[ISl . [ 21 1 . For electrons, some theoretical studies argue 
that the SOTR coupling parameter decreases with in¬ 
creasing wire confinement |22| , however later experimen¬ 
tal measurements [ 2 ^ for widths ranging from 1.18 ^m 
down to 210 nm suggest the opposite. For hh, one fol¬ 
low that the factor associated to the linear term behaves 
as theoretically predicted for QID electrons, while the 
one associated to the cubic term do not depend on the 
wire width 0. Moreover, recent magneto-transport 
measurements in Ge/Si core/shell nanowires 0 show 
that SOTR strength of holes behaves anomalously, as in 
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the Q2D case. We get motivated about the undoubted 
existence of competitor high-order terms on the quasi¬ 
momentum within the SOI-R model for QlD-hh systems, 
and the lack of a formal treatment for the Q2D-to-QlD 
transition pursuing an appropriate QlD-hh SOI-R model 
that fulfill the issues presented above. 

In this Letter, we revisit the SOI-R modelling for Q2D- 
hh systems and we manage to improve the development 
by Winkler et al. 0 We show that our representation of 
the SOI-R coupling parameter, support its subtle anoma¬ 
lous behavior, when one let the density grow. We have 
deduced this quantity explicitly, compare it with avail¬ 
able experimental data as well as with other numerical 
simulations, and stress the remarkably better agreement 
we have found. A similar procedure allow us to derive a 
model Hamiltonian for SOI-R in QlD-hh systems. 

It is widely accepted that SOI-R for the valence band 
(VB) is a second order effect. We are interested in quan¬ 
tum systems made of direct semiconductors with zinc- 
blende structure (point group T^), in particular those of 
III-V and II-VI binary compounds. Consequently, we 
start from the k-p (8 x 8) Pidgeon-Brown (PB) model in 
the spherical approximation |24l l25l| , which involves VB 
and the conduction band (CB). The inclusion of CB is 
also supported by results [uj which suggest that SOI- 
R is strongest in narrow-gap semiconductors. Then, by 
considering proper confining potential, the Q2D or QID 
are focused. 

To accurately describe a Q2D system of confined 
holes at the heterostructure interface, the Poisson and 
Schrodinger equations have to be self-consistently solved. 
In the framework of our scheme, we model the het¬ 
erostructure potential profile as a triangular well, thus 
V{z) = eFz, where F = eNs/{2eeo) [l^. The present 
approximation had been successfully applied in the study 
of spectral properties of electrons [iSi as well as for hole 
0. 

Focusing on extracting from the Schrodinger eight- 
component system [ 2 ^ 


HpB + {V{z) - f )I8 


^{r) = Os, 


( 1 ) 


an effective Hamiltonian for descri bing the hh states, we 
use the partition Lowdin scheme 11 .0 up to second 
order of perturbation 


rr(0) _ ^0 

^mm' 

rr(l) ^ tt ! 

^mm' 
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Em ~ El 


H. 


( 2 ) 




Tjf 


( 2 ) 

(3) 


in Ref. [ 2 ^ — which couple the heavy holes with dif¬ 
ferent total angular momentum projection. Taking into 
account that Fl^^, and F[^^, are both null 0, after 
some simple algebra in the term (|3]) we finally obtain 


Hhh-\,hhi = 

phh _ JZlh 
-^n 

E'J:^ - E^^ 


(5) 


and 


Efih^,hh^ — 

f^hh _ phh 
-^n 

+ 


]^hh _ J^ss 


( 6 ) 


In the equations ([5]) and (jS]), L and M are terms of HpB 
given by 


mo dz 


M = 


2mo ~ 


(7) 


where k± = E iky, h is the Planck constant, mg 
is the bare electron mass and 7 ^ = ( 271 " -|- 37|')/5 — 
Ep/{6Eg), with 7(^2 3 the semi-empirical Liittinger pa¬ 
rameters, Ep is related to the Kane parameter Pk = 
—i{h/mo){S\px\X) through Ep = {2mo/h'^)P^, and Eg 
is the energy gap. Likewise, 




hh,lh,ss / \ _ 


(z) = 


A 




0 

hh,lh,ss 


Cn 


( 8 ) 

are the normalized solutions of the VB transversal states, 
for the diagonal part of iTpe with a triangular quan¬ 


tum well 171, with z[ 


0 

hh^lh^ss 


{^'^hh,lh,ss^ 


E/tF)K 

'"'hhjihlss effective mass parallel to the growth di¬ 
rection z, for the heavy-, light- and spin-splitted split- 
off holes: mo/( 7 i - 27 s), mo/( 7 i -I- 2^s) and mo/( 7 i), 
respectively. We represent by A one of the Airy func¬ 
tions, while B —the other one— diverges at 00 , being 
Cn = 2.338,4.088,5.521,... its zeros. Furthermore, the 
transversal energy levels are given by 


with and being the diagonal and off-diagonal 
component of iTps, respectively. 

We are interested in the elements F[hht,hhi and 
Hhhi,hhi^ —Hsfi and iLg.a in JTpb for the basis chosen 


j^hh,lh _ 

-C/„ — Cn 


JtBL + (9) 


with Aso as the split-off energy. 
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Finally, the terms dS]) and ([6l) lead us to recover the 
ex pec ted form of the Rashba model for the Q2D case 


j?r-2D = /3(^+fc-+<^-fc+), (10) 


and to derive a new expression for the SOI-R parameter 

ml ^ ^ 

where (t± = \{crx i icTy) and cr^ y the Pauli matrixes, 

while and are Airy function integrals of the 

form 

JMM ^ . ( 12 ) 


It is worthwhile to underline that the non-perturbated 
term (diagonal part) includes the electric field and we 
were able to solve it exactly, which is a differentiating 
advantage of the present approach respect to earlier mod¬ 
els [13 , which considered the electric field perturbatively 
and only in the third-order terms can be approximately 
obtained the none-zero off-diagonal matrix elements. We 
had disregarded the third-order perturbation term due 
to its contribution of about 0.2 percent. 

The deduced expressions (0), ® and (EH) represent 
a reliable starting-point platform to quote SOTR cou¬ 
pling parameter (fTTIl for specific semiconducting Q2D- 
hh systems, playing for hh the same role as expression 
(13) reported before for electrons 0. The completeness 
of these self-complementary exact formulae is unprece¬ 
dented for holes, as far as we know. 

Following expression (EH), the existence of an inverse /3 
dependence on energy-levels ® gaps for the VB states is 
observed. This suggests that VB levels separation, which 
is in principle an electric-field F strength-assisted raising 
process and/or a charge-carrier concentration Ns growing 
function, governs the charge carriers confinement. Fur¬ 
thermore, our result EH), is consistent with the anoma¬ 
lous Rashba spin-splitting (iHEHl . 

Fig. [T]displays with a func¬ 

tion of Ns (solid line) following EH- Besides, pursuing a 
validation of our modelling approach, a comparison with 
some relevant analytical 121, experimental 10, U , 28 - 
and numerical 0,0 contributions, had been included. 
In the case of experimental measurements the samples 
are of the type: Ala;Gai_xAs/GaAs, for several values of 
molar composition x and acceptor doping. Before contin¬ 
uing, let us briefly mention some specific details, namely: 
(i) The expression (fTTIl is “autonomous” upon events out¬ 
side the region where the Q2D gas is built and depend 
solely on GaAs layer’s parameters, {ii) In the case of 
Refs. 28-33, reports of carriers’ concentration measure¬ 
ments with different spin polarizations (Ab|_,V_) were 



FIG. 1. (Color online) Q2D conpling parameter /3i for the 
first hh subband, taken after EH, is shown (solid line), and 
compared with previous analytical flH (d ashed line), experi- 
mental (symbols: ■ Ref. [T^, A Refs. (^. 1^. O Ref. , <) 
Ref. [iini ]- V Ref. [dll ], and ♦ Ref. ^3) s-nd numerical jlfll. lll| 
results (dotted line) based on a self-consistent method. 


taken, and had been considered an expression that relate 
them with /3i 0,0 




2 Ns{N+ - Nr) + ^Ns{N+ + Nr) 
6Nl + 2AN? 


(13) 


being = y'Ng ± ANg . Confronting our model EH 
to the above mentioned experimental data, general rea¬ 
sonable coincidence has been found. Notice the re¬ 
markably good agreement within the low density range 
[I X 10^° — 5 X cm~^, when our simulation (solid line) 

crosses almost all experimental-point uncertainty zones, 
following as well the overall experimental-data trend. For 
larger density section [20x I0^°—40x 10^°] cm~^, we have 
achieved a better match for /3i concerning its measured 
quantity, meanwhile a slight mismatch in tendency have 
been detected. Finally, we turn to previous numerical 
simulations 0 - 0 . The case of Refs. 00 , repre¬ 
sents a local very precise fitting (dotted line) as can be 
seen in the region [2 x 10^° — 3 x I0^°] cm~^. On the con¬ 
trary, in the range [3 x 10^° — 4 x 10^°] cm“^, that nu¬ 
merical prediction noticeably departs from experiment. 
Globally speaking, we think that this calculation 00 
do not fully reproduce the overall tendency of the exper¬ 
imental points, in the range [1 x 10^° — 5 x 10^°] cm“^, 
which contrast with the finest adjustment exhibited by 
the solid curve EH in the same range. Worthwhile to un¬ 
derline that analytically quoted results (dashed line) of 
Ref. 01, are clearly far both from experimental points 
and from the above discussed numerical modelling (in¬ 
cluding ours). However, the general trend is some how 
preserved, specially in the high density region. 

FigEl plots Pi within typical density range [10^° — 
10 ^^ cm~^] and exercised several semiconductors with 
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Ns xlO“ [cm-2] 


FIG. 2. (Color online) SOI-R conpling parameter Pi as a 
function of Ns, for the lower hh sub-band in different semi¬ 
conductors. 


narrow: InSb (0.235el^, dashed-double/dotted line), 
InAs (0.417el^, dashed-dotted line); middle: GaSb 
(0.812 el^, dotted line); and wide: InP (1.4236 eP, 
dashed line), GaAs (1.519 eP, solid line) energy gaps, re¬ 
spectively. One can easily observe all analyzed materials 
satisfy the predicted anomalous trend. Nonetheless, the 
expected llj increasing ordering of Pi with the energy 
gaps is unsatisfied for the InP in the considered density 
range. 


Lets turn now to the problem of a QlD-hh systems, 
which can be lithographically achieved by placing repul¬ 
sive electrodes on top of a heterostructure. Following the 
Ghesi et al. [isj] approximation, by substituting ~ Q, 
~ and ~ 0 in the Hamiltonian (fTOl) . for a 

channel with lateral extent aligned with the y axis, 
we get 


TThh 

^R-ID 


Stt^ 


PCTyky + PCTyky, 


(14) 


which is the expected dependence we want to obtain 
from our modelling. We simulate the additional con¬ 
striction by an harmonic potential ^nioui^x^. The same 
perturbative treatment described above, starting from 
the Pidgeon-Brown model in the spherical approximation 
and carried out for the Q2D case, enable us to derive an 
effective SOI-R Hamiltonian for QlD-hh systems 


n = -a<jyky + a'a-ykl. (15) 




Importantly, a further cornpetitor cubic contribution 
arises as predicted before As a bonus, an appeal¬ 

ing concordance with foretold signs is nicely reproduced 
both for fcy-linear and A:y-cubic terms. 


The SOI-R coupling parameters in (1151) are given by 
_ _3^ ( 

2 I Tphh TPlh -r^hh TPSS I ’ 

^0 V ^nn' - ^nn' ^nn' “ / 

(16) 

< \ - E±, - A-, ) ’ 

(17) 


where the subindex nn' stands for the transversal states, 
with eigenenergies 


Enn’ = ^V7l ± Is {n' + ^)+Cr, 


{heF)i 


Enn' = {n' + ^) + C, 


jheF) 

' ^2m* 


,y^^*hh,ih ( 18 ) 


■A 


so, 


where 71 = yf - Fp/{ 3 Fg). 

The coefficients involved in (|16l) - (fT7l) are those in (fT^ 
—as in the Q2D case— along with 


jhl/hs 


/ hs 




2 


(19) 

( 20 ) 


where 


lhh,ih = \ -V 71 ± 7s 

V rriQUJ 



( 21 ) 


In additions to the Airy functions ([ 8 ]), transversal 
eigenstates {x, z) = com¬ 

prise the properly normalized Hermite polynomials. 



A similar discussion as the one proposed for the Q2D- 
hh coupling parameter /3 (ED, is valid here. Both QDl- 
hh coupling parameters (IISD-(ED decrease with an in¬ 
creasing external or intrinsic electric field, as shown in 
Fig. [3l This trend is in qualitative accordance with 
the SOI-R strength measured in holes Ge/Si core/shell 
nanowires 0 , depicted in the inset of Fig. [3l We un¬ 
derline that we had even recovered the same order of 
magnitude. 

We observe a non-trivial dependence on the charac¬ 
teristic confinement length, = l^h as in ( 12111 . w hich 
departs from the one in the approximation (fT^ jl 8 l | . In¬ 
deed, at typical densities (W > 10 ^^ cm“^) whose chan¬ 
nel extent along 2 : direction is ~ 20 nm, a and oP 
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Ns xlO^® [cm"^] 



FIG. 3. (Color online) QID coupling parameters Qn (solid 
line) and a'n (dashed line) for GaAs, taken after (1161) and 
(I17II respectively. = Ihh = 11.8 nm. Inset: SOI-R 

strength (-•-) measured in 10 nm width holes Ge/Si core/shell 
nanowires from Ref. [l^ . 


FIG. 4. (Color online) QID coupling parameters Qn (solid 
line) and Qn (dashed line) for GaAs, taken after (1161) and 
(HU) respectively, for typical {Ns > 10^^ cm panel (a)) and 
low {Ns = 10^° cm“^, panel (b)) heavy hole densities. 


decrease whenever the QlD-hh system become wider, 
as shown in Fig. [4] (a). However, a asymptotically 
tends to zero, while a' promptly turns into a plateau 
at the vicinity of /3 d- This is intuitively clear and 
straightforward from (HH). In short, our modelling of 
the QlD-hh systems recovers the Q2D-hh case, in the 
limit when grows. On the other hand, at low densi¬ 
ties {Ns = 10^° cm'^jL^ Ri 20 nm), we found a critical 
confinement length ~ 11 nm where both coupling pa¬ 
rameters exhibit an essential asymptotic discontinuity, as 
shown in Fig. [l](b), which resembles a Fano-like profile. 

Owing to a deeper analysis, it seems that this behavior 
is related to the hole effective mass anisotropy. Indeed, 
in the spherical approximation, hh effective mass perpen¬ 
dicular to the growth direction = mo/ (71 -b 7 s)) 

is lighter than Ih effective mass (m;*^j_ = mg/ (71 — 7 s))- 
Note that the behavior of a and a' is governed by the 
denominators and in ([161) and 

dni). The latter hardly vanishes due to the split-off gap 
Aso- On the other hand, for the first hh and Ih subbands 
(n = 1 and n' = 0), following (fT51) we have 



where A = (ji -b 7 s - (71 + 7 s) ^ (71 - 7 s) ) and A = 
^(71 -b 27 s) 3 — — 27 s) Both terms in (1^ are 


positive definite, thereof there is a value of the charac¬ 
teristic confinement length 


Ll = 


A 

CiA 


2 


(27) 


that cancels AS^K For Lu, = , a and a' diverge, 

if > Ll are both positive definite while for L^j < 
Ll are both negative definite. Despite Shubnikov-de 
Haas experiments have been performed for hundred-of- 
nanometers wide quantum wires 23|, the beating pat¬ 
tern become considerably weaker for narrow enough sys¬ 
tems, where the discontinuity take place (L^, k, 10 nm 
for Lz ~ 40 nm). On the other hand, discontinuities may 
be detected in magnetic focusing experiments [mill for 
QPC patterned in [001]-grown Q2D hole gases, by ana¬ 
lyzing the shift of signal peaks as a function of SOTR 
strength. As we have shown in FiglU the last quan¬ 
tity can be tuned by Lw Therefore, for different QPC 
widths the signal-peak shift mentioned above, may be 
tailored as well. We foretell an interplay in the broaden¬ 
ing of the signal-peak separation, whenever some critical 
width (alike ) is reached, leading then to the oppo¬ 
site tendency for the peaks evolution having trespassed 
some typical confinement width. The last could sup¬ 
port the very existence of such ranges of discontinuity. 
Nonetheless, in QID systems aligned along the growth 
direction z, this discontinuity no longer arises, because 
of the effective mass isotropy {A£^^ never vanishes) in 
the transversal plane {xy), as may be inferred from stud¬ 
ies of s pec tral properties of holes in QID systems under 


SOI-R m 13J,13£ 


In summary, multiband k-p (8 x 8) standard Hamilto¬ 
nian MM , is a reliable candidate to address the SOI- 
R coupling modelling for heavy holes embedded in both 
Q2D and QID systems. While we recover the expected 
form and phenomenology of the effective Hamiltonians 
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181, our accurate SOI-R coupling parameters deriva¬ 
tion is confirmed by quantitative and qualitative compar¬ 
ison with experimental measurements. According to our 
modelling, further estimation procedures of these quan¬ 
tities in specialized semiconductor alloys are no longer 
needed. As far as Q2D-ft.ft. is concerned, we found to be 
in a remarkable better agreement with respect to previ¬ 
ous theoretical simulations. SOI-R coupling parameter 
anomalous behavior is straightforward understood from 
our closed analytical expression. We had formally de¬ 
rived an effective SOI-R Hamiltonian whose cubic term 
is competitor to the linear one. The coupling 

parameters also behaves anomalously, and follow recent 
experimental observations. In this case, we identify two 
different trends in the dependence with the characteris¬ 
tic confinement length of the nanowire. While at typ¬ 
ical densities both coupling parameters decrease as the 
system become wider, tending to the Q2D case, in the 
low density regime they exhibit an essential asymptotic 
discontinuity, that may be produced by permutation of 
hh — Ih effective mass roles. 
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